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S ■ Abstract 



In the present work the evolution of a coherent field structure of the Sine-Gordon 
equation under quantum fluctuations is studied. The basic equations are derived 



Q I from the coherent state approximation to the functional Schrodinger equation for the 

' field. These equations are solved asymptotically and numerically for three physical 

situations. The first is the study of the nonlinear mechanism responsible for the 
^ . quantum stability of the soliton in the presence of low momentum fluctuations. The 

second considers the scattering of a wave by the Soliton. Finally the third problem 
considered is the collision of Solitons and the stability of a breather. It is shown 
that the complete integrability of the Sine-Gordon equation precludes fusion and 
splitting processes in this simplified model. The approximate results obtained are 
non-perturbative in nature, and are valid for the full nonlinear interaction in the 
limit of low momentum fluctuations. It is also found that these approximate results 
are in good agreement with full numerical solutions of the governing equations. This 
suggests that a similar approach could be used for the baby Skyrme model, which 
is not completely integrable. In this case the higher space dimensionality and the 
internal degrees of freedom which prevent the integrability will be responsable for 
fusion and splitting processes. This work provides a starting point in the numerical 
solution of the full quantum problem of the interaction of the field with a fluctuation. 
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I. INTRODUCTION 



In the past few years it has been shown that the Skyrme model [|I| is related to the low energy 
limit of QCD p[. This fact, together with progress in the approximate and numerical solution of 
strongly nonlinear equations, has renewed interest in a detailed study of the quantum mechanics 
of the Skyrme model. These studies are directed along two main lines. The first simplifies the 
model at the classical level to the so-called Baby Skyrmion model, which in turn is related to the 
Sine-Gordon equation lll,^,^. In this simplified classical model extensive numerical studies have 
led to an increased understanding of the stability, scattering and interaction of coherent structures 
with radiation [0,^. On the other hand, the second line has focused on the quantum effects of 
the full Skyrme model P,^. In these studies the quantization was either obtained by linearization 
around a field configuration or by a finite dimensional approximation to the full problem ^ . 

Here we take a complementary approach. We study the 1+1 Sine-Gordon equation keeping all 
the degrees of freedom and quantize along the lines of [[7|-pil[|. This leads to a field equation which 
is strongly coupled with the equations for the fluctuations. In [jTO[ the formalism for the functional 
coherent state approximation was fully developed and the possible advantages and disadvantages 
of various approaches designed for numerical purposes were discussed in detail. In particular, 
the closed-time path method with the Hartree factorization (see e.g. [|lT|) was applied in |l^ 



to the static Sine-Gordon field, and the approximate results obtained for phase transitions and 
stability were found to compare favorably with known exact results. In this same work static 
results were also obtained for more realistic field equations. Our present work differs from the 
above cited papers in that we choose to approximate the Green function (the variance kernel 
for the Gaussian ansatz) by a suitable parametrized trial function. This choice leads to a great 
simplification of the problem for the case of low momentum fluctuations. Thus, for example, the 
infinite system of partial differential equations of become ordinary differential equations and, 
eventhough the field and the fluctuations in our formalism are strongly coupled, we are able to use 
some of the solutions of the classical Sine-Gordon equation (since in 1 + 1 dimensions it is completely 
integrable) to construct approximate solutions to the quantum problem, including the effect of the 
radiation generated by the quantum fluctuations. We then solve the equations numerically and 
these numerical solutions are compared with asymptotic solutions. Note also that the approach 
used in the present work is completely different to that of |[T^ , pr3[| where the wave functionals are 
constructed using Gaussian approximations to the functional Schrodinger equation for the Sine- 
Gordon field. However, the particles are considered as higher excited states (in function space) of 
the linearized field equations. In our treatment the field equations are nonlinear and dynamic and 
different particles are represented by different nonlinear field configurations, not by higher order 
Hermite functionals as in |]12|,|13|. 

Finally, it must be stressed that the approximate analytic results obtained here are valid in a 
strongly nonlinear regime and, in principle, do not depend on the 1 + 1 nature of the model and 
could be used to study problems related to more realistic simplifications of the Skyrme model. 

The paper is organized as follows. In Section II the detailed formulation of the quantum problem 
is stated, with the free parameters adjusted to mimic mesons and baryons. Section HI is devoted to 
the study of the coherent state approximation and the derivation of the quantum equations for the 
field for low momentum fluctuations. In Section IV three problems are considered. The first is the 
nonlinear stability of a single Soliton under the influence of quantum fluctuations. This stability 
is studied both numerically and asymptotically. In particular the asymptotic solution includes 
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the damping effect of the radiation shed by the Sohton due to the ffuctuations. This asymptotic 
solution explains in detail the mechanism for the nonlinear stability of the Soliton. The second 
problem studied is the scattering of a meson (wave) by a static Soliton. The numerical solution 
for this problem shows that in this process the Soliton is also stable. Finally the third problem 
studied concerns the collision of Solitons and the quantum evolution of a bound state (Soliton and 
an anti-Soliton). It is shown that the complete integrability of the Sine-Gordon simplification of 
the full Skyrme model to just one internal degree of freedom for the field precludes the processes 
of fusion and splitting. The processes of fusion and splitting are produced by the influence of all 
the internal degrees of freedom. 



II. FORMULATION OF THE PROBLEM 



For the basic structure we take the baby Skyrme model, which is a reduction of the full model 
with only two fields present This model, in turn, reduces to the Sine-Gordon equation which 

is known to be completely integrable H]. In these variables the Hamiltonian takes the form 



H = mc I 



-TT 



1 / V 

2 \dx 



+ I ^(1 — COSip) 



dx, 



(2.1) 



where ip is an angle variable whose momentum vr is given by 



TT 



Idip 
c~dt' 



(2.2) 



m is the mass of the particle and / is a typical particle size. Using dimensionless variables x = x/l 
and t = ct/l, we obtain, after dropping the tildes. 



H 



mc 



-IT 

2 



1 f^A 

- I — I + (1 - COSV9) 



dx. 



2 \ dx^ 

with TT = dif/dt. To mimic a baryon by means of the Sine-Gordon soliton, we take I ~ 10"^^ 
and m ~ 10"^'' Kg, which gives an internal time of 10~^^sec. 

The equation of motion derived from the Hamiltonian variational principle 



(2.3) 



cm 



fti 
6 / (tt^ 



H) dt 



is the Sine-Gordon equation 



+ sm.ip> = 0. 



The one dimensional Skyrmion is the soliton solution 



-4 arctan 



exp (^—{x — vt) I Vl — 



(2.4) 



(2.5) 



(2.6) 



of |2.5| , which represents a localized deformation at x = vt. In this context v must be taken to satisfy 
V <C 1, since the Skyrme model is only consistent for small energies. The linear travelling periodic 
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wave solutions of the Sine-Gordon equation are interpreted as pions. Notice that the nonhnear 
model contains linear waves which describe bosons and nonlinear structures which describe fermions. 
In his original work Skyrme suggested that the fermionic part of the Lagrangian is needed only 
just to count the number of localized states of finite amplitude of the bosonic field [0]. In this 
interpretation, fermions are just a point approximation to nonlinear localized bose fields. It has also 
been suggested on rather general grounds that the canonical quantization of fields (as bosons) gives 
Fermi-Dirac type statistics for the kinks |TB[. In the Sine-Gordon model the exclusion principle 



holds for kinks, since we know that for the general exact solution there is no solution with two 
identical solitons \Ti\. 



In this article, we shall consider a canonical quantization using the functional Schrodinger 
picture. It is to be noted that all quantizations for the Skyrme model cited in the Introduction 
make the same assumption. However, in this work we differ from previous treatments in that we 
shall keep infinitely many degrees of freedom in the classical field and reduce the dimensionality 
of the space of fluctuations. The final result will be shown to be a system which consists of a 
partial differential equation (similar to the Sine-Gordon equation) for the field coupled to a system 
of nonlinear ordinary differential equations which control the fluctuations. 

In the field configuration representation the Schrodinger equation takes the form 

a'JM = Hm<p), (2.7) 

where '${(p{x),t) := (y9(x)|\l/(t)) is the amplitude for finding the field system characterized by the 
state vector |^E'(t)) in the field configuration v^(x) at time t. In this configuration representation 
the scalar product of two state vectors is given by the functional integration 



and the field operators are represented by functional kernels. Thus the field operator 0(x) is 
represented by {ip{x)\(j){x)\'^ (t)) = ip{x)'^{{p{x),t) and, therefore, acts as a multiplication operator, 
while the action of the canonical momentum operator is given by 

{V{x)\n{x')\m) = -^nJ^^{v{x)Mt))■ (2.9) 



The Hamiltonian field operator H then becomes 
H = mc 



dx, (2.10) 



where there is no ambiguity in the ordering and the functions of operators are defined by their 
power series. 

The quantum mechanical problem for the field consists of solving the Schrodinger equation ( p.7|) 
for a given initial field configuration. Notice that the time in ( p.7|) has a scale set by mc^ /%, which 
is of the same order of magnitude (10~^^) as the time scale c/l for the field fluctuations. This is 
to be expected since the Skyrme equation was assumed to be consistent with the quantum scale of 
the particle. Thus we can take the same dimensionless time variable for either scale. 
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III. APPROXIMATE SOLUTIONS TO THE FUNCTIONAL EQUATION 



To solve the Schrodinger equation ( |2.7| ), we take a coherent state approximation [0-[T^ and 
study the evolution of its parameters. Following [ffHTOH we consider the functional 



which we extremize with the Gaussian trial functional 

t) = exp |i [ 7i{x, t) [^{x) — t)] dx 
exp 



(3.1) 



dxdy ['{){x) — 0(x,t)] 
y, t) - iS(x, y, t) [^{y) - ^{x)] \ , 



(3.2) 



where the kernels and S take into account the quantum fluctuations and (l){x,t) and TT{x,t) 
are the average field and average momentum of the Gaussian, respectively. 

Substituting the trial function ( ^.21 ) into the functional (^) and integrating over the (p{x) 
variable yields an averaged action [[7|-[T0[1. This action is given in terms of (j){x,t), Q and S. The 
potential is then also expanded around the average 0(x, t). It should be stressed that the consistency 
of this approximation depends on the smallness at all times of the variance Q, in the sense that 
the average energy of the fluctuations around the mean 0(x, t) is small compared with the energy 
of the average motion. Moreover, by keeping only quadratic terms resulting from the averaging 
around the mean of the nonlinear term, which amounts to making the assumption that the energy 
of the fluctuations is small compared to the energy of the mean, we arrive at an effective action of 
the form 

dip 1 fd(p 



(3.3) 




+ TiCl{x, X, t) — 2Sr2S(x, X, t) — cos{ip)Q{x, x, t) 

dxdt + 0{n^) 



+-^r^{x,y,t)\^=y- -il {x,x,t) 



2 dx^ 



where the notation fiS is taken to mean the kernel of the operator defined by the convolution of E 
with Q, that is 



QT.{x,y) 



Q{x, z)S(z, y) dz 



and SfiS is the kernel of the operator defined by the convolution of S, and S 



SfiS(x, y) 



-L J-L 



S(x, z)Q{z, u)T,{u, y) dudz. 



(3.4) 



(3.5) 



In the following we will take the limit L ^ oo at different stages in the calculation of the effective 
action. The effective action can be computed once we choose an appropriate parameterization 
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for the variance and the phase S. It must be noted that thanks to the simple form of the 
potential, the Gaussian integral can be evaluated exactly. However, since we are interested in 
small fluctuations we stop at the quadratic level. We will come back to this point when comparing 



our results with those in 10 



Now, since the field is homogeneous, we can take 

1 

fl(x,y,t) = — / exp[ik(x — y)]il{k,t) dk (3.6) 

T,{x,y,t) = — / exp[ik{x — y)]T,{k,t) dk, (3.7) 
2n J-oo 

where 

The parameters a{t) and b{t), which control the spreading of the fluctuations, are to be deter- 
mined after the effective action is varied. The choice of the trial function is then guided by the 
simplicity of the resulting expressions. Notice, however, that since the results obtained depend 
only on the spreading, the same qualitative behaviour will be obtained for other forms of the trial 
function. 

Since the approximate solution ( p.2|) involves the kernel the proposed expression is con- 
vergent provided that the momentum of the fluctuations involved in the integration is low. This 
assumption is in agreement with the fact that the basic Skyrme model is consistent at low mo- 
mentum and the assumed homogeneity of the fluctuation. This is taken into account by taking for 
Q~^{x, y, t) the cut-off kernel 

n-\x, y, t) = e''^^"-^) f k^ + a") dk. (3.10) 

2-nilQ J-K ^ ' 

Then the non-constant contribution of fi~^(a;, x,t) is 

In a similar manner we obtain 

1 52 , , I ^ k^ 



o i k 

nTr^^i^^y^'^)\y=^= ~^^o / TT- — i^dk 
2 ox^ K^oo An J-K k^ + 

ftoK 1 dk 



lim 



+ -a{t)Qo J_ 



2tt 47r ' ' 7-00 fc2 + 1 



(3.12) 



The first term in this expression is infinite, but a constant. Thus it does not contribute to the 
equations of motion. We therefore take just the second term in the effective Lagrangian. Finally 
the parameter Qq measures the size of the fluctuations. In principle, other choices of the parameters 
may lead to different regularizations. However, as discussed below, the basic qualitative picture 
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described in this work is not changed by these alternative regularizations, provided the momentum 
is low. For the case of higher momentum, procedures similar to the ones discussed in [^,10| can be 
used. 

With the above assumptions, the effective Lagrangian takes the form 



^ \ d'-P 
J-L I dt 



1 / V 

2 \dx 



+ 



2Txa 



C2 cos 1^ 



where the constants Ci and C2 are given by 



a + -\lQC2a I > axat, 



(3.13) 



Cl 



dk 



IT 



[l + k 



2\3 



C2 



dk 



TT. 



(3.14) 



It is now convenient to change variables and define q = ^ and b = p in order to obtain the 
Lagrangian ( p.l3|) in the form 



2 



+ 



TC- 



'dt 



-TC 



1 / V 

2 V9x 



il 



2 Cl 
2tx 



L 
— ( 



Txq 3 ) cos if] 



dx > (it. 



(3.15) 



The equations of motion are obtained by varying the Lagrangian (|3.15|) with respect to the param- 
eters TT, yj, p and q. These variational equations will consist of a partial differential equation for 
the field coupled to ordinary differential equations for the fluctuations. For the field the variational 
equation is 



The equations for p and q are derived from the Lagrangian 



2ciLfio 



[qp - 3H{p, q)] dt, 



3tt Jo 

where the Hamiltonian for the fluctuations is given by 

1 



H{p, q) = p^q + q 



,1/3. 



2L J-L 



L 1 /■ 

cosy? dx — -q^^ 



K 



-2/3 



(3.16) 



(3.17) 



(3.18) 



The equations of motion for p and q are then given by Hamilton's equations as 



■ dH ^ „5/3 
q = — = Qpq 
op 



p 



OH 
dq 



5p q 



2-8/3 



+ 



-5/3 _ J_n-^/3 _ I 
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q 



'2/3 

2ir 



cos (p dx, 



(3.19) 
(3.20) 



where 



/5 



3K 



Since K is assumed to be small, we will take (3 
section. It must be noted that for initial conditions which have 
X ^ oo, or vice versa, 



(3.21) 

2.5 in the numerical calculations of the next 
27r as x ^ — oo and — > as 



A 



lim — 



COS (f dx = 1. 



(3.22) 



It is therefore apparent that the Equations (|3.2CI|) for the fluctuations decouple from the Equation 
( CT ) for the field. 

It is interesting at this point to compare our equations (|3.16| ), (|3.19|) and (|3.2CI| ) with the cor- 
responding equations (4.5) of Ref. [10]. Observe that in that reference the equation for the field 
takes the form 



Taking [3=1 and assuming G{x, x, t) ^ 1, 
form of small fluctuations, we have 

„-iG(x,x,t) 



^^(^'^'*)sin/3v? = 0. (3.23) 
which amounts to choosing the initial conditions in the 



1 



^G{x,x,t). 



(3.24) 



Clearly, substituting this last expression into ( p.23| ) recovers our equation ( p.l6| ), so the respective 
field equations agree for the above mentioned initial conditions and value of (3. Now, as for the 
second equation (4.5) in 



TDI 



note that this is an infinite system of partial differential equations 
for the operator valued function ^, while in our formulation, because of the assumption of spatial 
homogeneity and the functional form chosen for the trial Green function, the system simplifies to 
a problem of ordinary differential equations. Another difference between our approach and that 
followed in ||10|| is that the equations proposed there in order to arrive at numerical solutions are 
integrodifferential equations, as opposed to the simple variational approximation we propose for 
obtaining solutions. It must be also remarked that our assumed homogeneity of the Green function 
is consistent with the low momentum limit we have chosen. In this low momentum limit, the 
fluctuations do not resolve the fine scale of the field and, to leading order, the configuration is a 
homogeneous background for the fluctuation. To conclude this section we also consider important to 
stress the fact that in previous published work the interest has been on static solutions allowing for 
arbitrary momentum of arbitrarily large fluctuations. This leads to a different renormalized version 
of the gap equation [0,H,|1O| and, for large coupling /?, to a loss of stability (phase transition). Our 
approximation does not capture this region, since we have assumed from the onset small fluctuations 
and small momentum. However it must be noted that our procedure could be extended to handle 
large momenta by choosing different trial functions for G and S, similar to the ones used in PJ. Also, 
due to the special form of the potential in the equations, the Gaussian integral may be evaluated 
to a better degree of approximation, thus allowing to include fluctuations of a larger amplitude. 
This program is currently under investigation and will be reported subsequently. 

In the following section we undertake a detailed study of the dyamics described by the equations 
( CT ). O and (|Oq ). 
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IV. SOLUTIONS 



The system of the Sine- Gordon equation ( p.l6|) and the equations (|3.2CI| ) and ( p.l9| ) for q and p 
describe the couphng of the field to the fluctuations and the corresponding feedback. Note that the 
fluctuations have been assumed to be small. However, they are allowed to feed back onto the basic 
field configuration. We shall now use these equations to describe in a nonperturbative manner the 
nonlinear evolution of some special field configurations. 



A. Quantum Stability of the Single Soliton 



We begin by studying the stability of the Soliton solution (p.6|) under a class of initial values 
for p and q. Note that a small value of q represents a small variance. Stability is then assured in 
the model by the fact that q remains small and that the field maintains its identity as a localized 
structure. 

Numerical integrations of the Sine-Gordon equation (|3.16| ) and the Equations ( |3.20|) and (|3.19|) 
for p and q have been performed for a wide range of initial conditions and typical behaviors are 
shown in Figures |l(a)| and |l(b)| . In Figure |l(a)| a numerical solution for $ = yj^: for the Skyrmion 
is shown and in Figure |l(b)| the behavior of a, the maximum of $, is shown. It can be seen that 
the fluctuations of the Skyrmion produce radiation, but that the field eventually stabilizes. This 
can be clearly seen from the maximum behavior shown in Figure |l(b)| . The stabilization onto a 
modulated small oscillation of the Skyrmion amplitude can be clearly seen. These results exhibit 
the strong stability of the Skyrmion with respect to fluctuations. It is possible to understand this 
behavior by making use of the modulation theory given in [|T^ and [|r^ by means of the following 
argument. If the scales for p and q are slow, we may take as an approximate solution 



Lp = Lp 



T ^0 1. 

1 03 

2 ■ 



which satisfies 



COS(p 



1 9 

2^ 



(4.1) 



(4.2) 



Hence as L — oo we have from ( p.22| ) that A = 1. With this, the equations for p and q are derivable 
from the Hamiltonian 



Hip, q) = 3 [p'q-^ + g^/^ _ ^ ^-1/3 + _^ 

Thus the orbits in the (p, q) plane are just the level lines of 

H{p,q)=E. 
The orbits of the (p, q) system are then given by 

5/3 



-2/3 



(4.3) 



(4.4) 



p = q 



E 



4U 







q'/' [E - V{q)] . 



(4.5) 
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The potential V{q) has a minimum which gives an oscillatory solution for p and q, so that the 
width of the Skyrmion and thus the amplitude of ^ = (f^ oscillate in time. The numerical solution 
shown in Figure 1 shows that the radiation, which is not taken into account in this approximation, 
stabilizes the oscillations onto a limit cycle. 

This strongly nonlinear mechanism accounts for the stability of the Skyrmion. In fact it is the 
feedback of the field on the fluctuations which produces the term g^/^ in V{q) and it is this term 
which stabilizes the motion. The potential V{q) has a maximum for small q. For energies E larger 
than this maximum, the fluctuations q increase and the field structure is destroyed. However the 
value of q for this to occur is too small to be consistent with the coherent state approximation. 
The model is therefore self-consistent and provides an explanation of how nonlinear interactions 
are responsible for the quantum stability of the field. 

The approximate solution above does not take into account the radiation produced by the 
oscillating Skyrmion and so this approximate solution will not give the baryon settling onto a limit 
cycle solution. To take account of the radiation the ideas of Smyth and Worthy |TB[ can be used. 
In this work the effect of shed dispersive radiation on the evolution of a single pulse for the Sine- 
Gordon equation was treated. To take account of the radiation we proceed as in [0], indicating 
only the main differences from this work. 

The Lagrangian density for the Sine-Gordon equation is 

' \'t-l^l-(l-%'^')il-cos^). (4.6) 



2^' 2 V 2 
To obtain an approximate solution of the Sine-Gordon equation, the trial function 

= -4arctane"^/"'(*\ (4.7) 

which is a soliton-like pulse with varying width w{t), is substituted into the averaged Lagrangian 



oo 



L dx 



oo 



^!£:_l_4(l-^l^,V3,,, (4.8) 

3 W W \ 2 



as in . In this approximation the Hamiltonian for p and q again does not change due to ( p.22|) 



The effect of the radiation shed by the evolving soliton is determined by finding an appropriate 
solution of the linearized Sine-Gordon equation . The effect of this radiation is then to modify 



the Euler-Lagrange equation for w{t). It is noted from the numerical solution of Figure P"(a]] that the 



radiation ip is of small amplitude compared with the soliton. Therefore following [|18| we consider 
the linearized Sine-Gordon equation 

for the radiation ip. This equation is solved together with appropriate source conditions at the 
pulse at X = 0. Since fio <^ 1, 



dt\ 2 
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It is then possible to obtain an expression for the radiation by making the adiabatic approximation 
that 1 — ^^q^/^ is constant to leading order. The effect of the radiation can then be found from the 
expression of |T^ by a suitable re-scaling. In this manner we find that the equations governing the 
evolution of the soliton, including the effect of radiation, are 



27r2 (Pw 



3w dt"^ ?>w'^ \ dt 



4 f 1 - ^gV3 



1 r 27r2 dw 
V\ [ SwVXt dt 

+ I Ji(v^(t-r))^^ 



dq 

'di 
dp 

'di 



3\'wt Jo 



dr 



(4.11) 



5p q 



2^-8/3 



+ 



-5/3 _ -2/3 , t^-ifi 



where 



A 



1 _ ^^0^/3 

2 ^ 



(4.12) 



These equations were integrated numerically. Comparisons between solutions of these equations 
and the full numerical solution of the Sine-Gordon equation for the amplitude a of $ = ip^ and 
q{t) for the fluctuations are shown in Figure 0. It can be seen that the amplitude agreement shown 
in Figure ^ is good considering the assumptions that were made to incorporate the radiation loss 
in the approximate equations. It can be seen that the approximate equations provide a suitable 
approximate solution for the full field behaviour using a finite dimensional approximation which 
includes radiation. It is noted that, since q{t) is periodic, the Sine-Gordon Equation ( p.l6| ) is subject 
to a parametric excitation. However the nonlinearity and radiation loss provide the necessary 
damping to enable a limit cycle to be achieved. 



B. The collision of a wave with a static soliton 

As a final example we consider the scattering of a wave packet representing a pion with mo- 
mentum k with a static soliton, representing a baryon originally at rest. The problem is solved by 
numerically integrating the Sine-Gordon equation ( |3.16| ) using the initial condition 

ip{x) = -4arctane"'' + f{x) (4.13) 

^ = 9{^)^ (4-14) 

where the functions / and g are given by 

f{x) = asmk{x + xo), \x + xo\ < 6 (4-15) 
g{x) = —ayW+lcosk{x + xq), \x + xo\ < 5. (4.16) 
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This initial condition represents an incoming meson with momentum k impinging on a nucleon 
located at x = 0. A numerical solution for the scattering of the pion wavepacket can be seen in 
Figures P(a)| to P(c)| . The initial condition (at t = 0) is shown by the solid line in Figure |3(a)| . 
In this Figure a reflected wave packet, a reorganized field configuration and a new packet shed by 
the baryon as a result of the interaction can be seen. In Figures |3(b)| and |3(c)| the complicated 
evolution of the baryon amplitude is displayed. This amplitude behavior is due to the interaction 
of the baryon with the packet. The scattering then involves a reorganization of the field, which is 
not taken into account when the particles are taken to be point particles. The description of the 
interaction of the baryon with radiation using a multi-phase solution of the Sine-Gordon equation 
is under investigation at present. 



C. Collision of two solitons in the presence of a fluctuation 



Since the classical field equation is completely integrable, solitons interact elastically and do 
not change configuration. The effect of quantum fluctuations on the collisions of solitons and this 
clean interaction will now be studied. 



Figure [4(a)| shows the collision of two Solitons with equal and opposite velocity. The initial 
condition used was 



27T - 4arctane-(^+^»-^*)/^'^ - 4 arctane-(^-^°+^*)/^'^ 



(4.17) 



as t — > — oo. Since there is no classical solution with twice the baryon number and zero velocity, the 
effect of the quantum fluctuations is just to slightly modify the classical interaction. The Solitons 
again settle down to a limit cycle for which the parametric resonance is balanced by the radiation 
damping. 

Figure |4(b)| shows the collision of a Soliton and an anti-soliton. The initial condition is 



-4 arctan 



a 



sinh 



It 



cosh ax 



(4.18) 



Again this interaction does not produce disintegration, just a modification of the classical interac- 
tion. 

Finally the susceptibility to disintegration of the breather-type configuration 



ip = —4 arctan 



a 



sin a/1 — t 



cosh ax a/1 — Q;2 



(4.19) 



with frequency Vl — is studied. From the numerical solution shown in Figure |4(c)| it can be 
seen that the breather is stable with respect to quantum fluctuations. 

The solutions displayed in Figure 4 show that the reduction of the Skyrme model to the Sine- 
Gordon equation is too severe for treating collisions. In order to obtain non-trivial collision and 
fusion processes, such as those possible for the nonlinear Schrodinger equation, reductions of the 
Skyrme model which retain more internal degrees of freedom must be derived. 
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V. CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 



We have formulated the quantum field problem for the Sine-Gordon equation which is related to 
the (dimensionally reduced) Skyrme model. Using the coherent state approximation for the solution 
of the functional Schrodinger equation, we obtain a solution of the partial differential equation for 
the (quantum corrected, semi-classical) field, which is coupled to ordinary differential equations for 
the fiuctuations. Other quantizations for nonlinear fields keep only finitely many degrees of freedom 
(minisuperspace approximation), which are then quantized in a canonical way. 

The first problem considered in the present work was the stability with respect to quantum 
fiuctuations of a Soliton. Both numerical and asymptotic solutions were considered. It was found 
that the nonlinear saturation of the field equation together with the loss of radiation balanced the 
parametric excitation of the fiuctuations. The fiuctuations in turn were controlled by the shape 
of the field. The good agreement found between numerical and asymptotic solutions suggests that 
finite dimensional approximations to the dynamics of the Skyrmion model, such as those used in 
1^, are also good approximations to the full dynamics of more complicated problems, such as those 
treated in 0. 

The scattering of a wave by a static soliton was also studied. The numerical results obtained 
show well defined waves and a Skyrmion after collision, which suggests the possibility of using 
multi-phase solutions, such as those of [^, to understand this scattering process. 

Finally several collision processes were studied. It was found that the reduced Skyrme model can- 
not account for the collision and fusion of baryons. Therefore the study of the fusion of Skyrmions 
into a toroidal configuration requires a uniform solution which interpolates between the torus and 
the individual Skyrmions. The possibility of using the solutions given in pO|j2T[| is currently under 
study. It must be noted that more sophisticated numerical formulations such as the ones proposed 
in [|10| must produce, in the limit of low momentum, solutions comparable to our results. 

To conclude, we note that the techniques described in this work can be applied to the study of 
low dimensional black holes. Indeed, an old observation that Sine-Gordon theory and 2-dimensional 
spaces of constant curvature are very closely related has recently found an interesting application 
to gravity in 1 -|- 1 dimensions. More precisely, Gegenberg and Kunstatter [^] have noticed that 
when a two dimensional Lorenzian metric is parameterized as 

ds^ = - sm^{u/2)dt'^ + cos^{u/2)dx^, (5.1) 

then the condition of constant curvature is equivalent to the condition that u satisfies the Euclidean 
Sine-Gordon equation. On the other hand, the so-called Jackiw-Teitelboim theory in two dimensions 

I = j <P{R- K)^/^dtdx. (5.2) 

has as solutions space-times of constant curvature R = A. Furthermore, the one-soliton solution 
of the Sine-Gordon equation has been found to represent (a patch of) a black hole solution of the 
(Jackiw-Teitelboim) theory [^. That a constant curvature space-time can be interpreted as a 



black hole is not unique to two dimensions. The 1-1-1 Jackiw-Teitelboim black hole can indeed be 
interpreted as a dimensionally reduced BTZ (non-rotating) black hole and many of its properties 
(including thermodynamics) have been studied ||23[ . 

To perform an analysis similar to the one presented in the present work for the Euclidean Sine- 
Gordon equation is cumbersome, since the equation is now elliptic and does not accept a well-posed 
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initial value formulation. However it is possible to work in the framework of a well-posed problem 
if one chooses a different parametrization for the two dimensional space-time as follows 

ds^ = - sin\i^{u/2)dt^ + cosh^{u/2)dx^ . (5.3) 

In this case, the constant curvature condition reduces to the Lorenzian Sinh-Gordon equation. It is 
then possible to analyze the quantum stability of a black hole solution using the functional methods 
presented in this article. This work will be reported elsewhere. 

As a final remark we point out that the quantum equations for a classical field obtained using the 
functional Schrodinger equation and the coherent state approximation will always have the same 
structure. Namely the classical equations for the field with renormalized (fluctuating) parameters 
and equations for the (parameters of the) fluctuations which are non-local in the flelds will always 
be obtained. 
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FIGURE CAPTIONS 



Figure 1. Stability of a single Soliton to fluctuations. Solution of Sine-Gordon equation ( p. 161 ) and 
equations ( p.20|) for p and q with Qq = 0.6 and (3 = 2.5. The initial conditions are q = 1.0 
and p = 0.0 and i; = in the soliton solution (|2.6| ) at t = 0. 

(a) Soliton. — : initial condition; : Soliton at t = 100. 

(b) Evolution of maximum a of $ = <^x- 

Figure 2. Stability of a single Soliton to fluctuations. Comparison between the full numerical solution 
of the Sine-Gordon equation ( 3.16 ) and equations ( 3.20|) for p and q and the approximate 
Equations ( [4.11| ). = 0.6, (3 = 2.5. The initial conditions are q = 1.0, p = 0.0 and f = at 
t = 0. Amplitude a of $ = (px- Full numerical solution: — ; approximate solution: . 



Figure 3. Scattering of a wave packet (pion) with a baryon. Solution of Sine-Gordon equation (|3.16|) 
and equations ( |3.20| ) for p and q with Qq = 0.6 and (3 = 2.5. The initial conditions are given 
by (^13D , (gmi) and a = 0.1, k = 1.0, 5 = 4.0 and Xq = 20 in ( ^151) and ( ^161) . Also q = 1.0 
and p = 0.0 at t = 0. 

(a) Solution at t = 50. 

(b) Evolution of maximum a of ^ = (px- 

(c) Detail of evolution of maximum a of ^ = (p^ from t = 15 to t = 40. 

Figure 4. Collisions of Solitons. Initial conditions have q = 1.0 and p = 0.0. Qq = 0.6. 

(a) Two Solitons. Initial condition ( [4.17| ) with xq = 15 and v = 0.2. Initial condition (t = 0): 
— ; solution at t = 150: . 

1.2. Initial condition 



(b) Soliton and an anti-Soliton. Initial condition ( [4.18 ) with a 
{t = —15): — ; solution at t = 15: . 

(c) Bound state of a Soliton and an anti-Soliton. Initial condition ( [4.19| ) with a 
Initial condition {t = —5): — ; solution at t = 45: . 



0.98. 
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